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Abstract. - We investigate particle production in heavy-ion collisions at RHIC energies as function 
of incident energy, and centrality in a three-sources Relativistic Diffusion Model. Pseudorapidity 
distributions of produced charged hadrons in Au + Au and Cu + Cu collisions at ^/snn = 19.6 
GeV, 62.4 GeV, 130 GeV and 200 GeV show an almost equilibrated midrapidity source that tends 
to increase in size towards higher incident energy, and more central collisions. It may indicate 
quark-gluon plasma formation prior to hadronization. 
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Introduction. — The precise calculation and prediction of transverse momentum and 
rapidity distributions of produced particles is of basic importance in relativistic heavy-ion 
physics. In this Letter we propose nonequilibrium-statistical methods [1] to investigate 
| analytically the gradual thermalization in rapidity space occuring in the course of particle 

_C ■ production at the highest available energies. The approach is tailored to identify the fraction 

of produced particles in local thermal equilibrium from their pseudorapidity distribution 
functions in heavy systems, with a focus on Cu + Cu and Au + Au. It may yield indirect 
evidence for the extent and energy dependence of a locally equilibrated parton plasma. 

There exist other theoretical approaches that allow to compute rapidity distribution func- 
tions for produced particles, albeit with less precision. Some of them are based on QCD, such 
as calculations within the framework of the Parton Saturation Model [2] . Ideal hydrodynam- 
ics is well developed in its applications to relativistic collisions, but more realistic dissipative 
hydrodynamic approaches are still in the early stage of theoretical development [3] . 

Thermal models are outstanding in their ability to correctly predict particle abundance 
ratios at midrapidity, or momentum integrated [4,5]. But since these approaches do not 
deal with nonequilibrium-statistical effects, one can not expect precise results for distribution 
functions whenever hadronic or partonic thermalization processes through multiple collisions 
on an event-by-event basis are important [5]. Within a thermal model, such effects could 
be simulated to some extent by different values of local temperature and chemical potential 
when investigating particle production at different rapidities. 

Hence, nonequilibrium statistics is the natural choice for a detailed description of the 
gradual approach to statistical equilibrium in relativistic collisions of heavy systems. Our 
Relativistic Diffusion Model (RDM) underlines the nonequilibrium-statistical features of 
high-energy heavy-ion collisions, but it also encompasses kinetic (thermal) equilibrium of 
the system for times that are sufficiently larger than the relaxation times of the relevant 
variables. 
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It is of particular interest in relativistic collisions of heavy systems to determine the 
fraction of produced particles that attains - or comes very close to - statistical equilibrium 
with respect to a specific macroscopic variable, such as rapidity. In the three-sources RDM, 
these are the particles produced in the midrapidity source. Hence we analyze Au + Au 
and Cu + Cu pseudorapidity distributions of produced particles at RHIC energies from 
^/snn— 19.6 - 200 GeV corresponding to beam rapidities 2/ma£c=3.04 - 5.36. We determine 
the transport coefficients and numbers of produced particles in the midrapidity source as 
functions of the incident energies, and centralities. 

Relativistic Diffusion Model. — In the Relativistic Diffusion Model, the rapidity 
distribution of produced particles at RHIC energies emerges from an incoherent superpo- 
sition of the beam-like components that are broadened in rapidity space through diffusion 
processes, and a near-equilibrium (thermal) component at midrapidity that may indicate 
local quark-gluon plasma (QGP) formation. 

The time evolution of the distribution functions is governed by a Fokker-Planck equation 
(FPE) in rapidity space [1,6-10] 
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-[R(y,t)r = --[j(y)[R(y,t)r\ + D y — [R(y,tW (1) 

with the rapidity y = 0.5 • \n((E + p)/(E — p)). The rapidity diffusion coefficient D y that 
contains the microscopic physics accounts for the broadening of the rapidity distributions. 
The drift J{y) determines the shift of the mean rapidities towards the central value, and 
linear and nonlinear forms have been discussed. 

Here we use \i = 1 (due to norm conservation) and v = 2 — q with the nonextensivity 
parameter [11] q = 1 corresponding to the standard FPE, and a linear drift function 

J(y) = (yeq-y)/T y (2) 

with the rapidity relaxation time t v , and the equilibrium value y eq of the rapidity [1,10]. 
This is the so-called Uhlenbeck-Ornstein [12] process, applied to the relativistic invariant 
rapidity for the three components Rk(y,t) (fc=l,2,3) of the distribution function in rapidity 
space [1,6,9, 13] 
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Since the equation is linear, a superposition of the distribution functions [1, 6] using the 
initial conditions Ri^(y, t = 0) = 5(y±y max ) with the absolute value of the beam rapidities 
Vmaxt an d Rd{y,t = 0) = 5(y — y eq ) yields the exact solution. In the solution, the mean 
values and variances are obtained analytically from the moments equations. The equilibrium 
value y eq is calculated at each centrality from energy- and momentum conservation among 
the participants. For symmetric systems it is y eq = independently of centrality, but its 
deviation from zero is important for precise calculations of rapidity distributions in case of 
asymmetric systems [10]. 

Pseudorapidity distributions. — If particle identification is not available, one has 
to convert the results to pseudorapidity, rj — — ln[tan(0/2)] with the scattering angle 9. The 
conversion from y— to 77— space of the rapidity density 

dN p dN dN 

-d^-E^ = J ^ {m)/{PT)) ^ (4) 
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Table 1: Produced charged hadrons as functions of centrality in Cu + Cu and Au + Au collisions at 
y^ijvjv = 19.6, 62.4, 130 and 200 GeV in the Relativistic Diffusion Model. The number of produced 
charged particles at each centrality is iV*£*, the percentage of charged particles produced in the 
thermalized source is n^. The ratio n n t/r y determines how fast the system of produced charged 
particles equilibrates in rapidity space. The widths of the peripheral sources are I\2, the width of 
the midrapidity source is r e q. The ^ j d.o.f . values with constraints are discussed in the text. 



System 


CTtot 


Nit 


n2(%) 


Tint 
T V 


ri, 2 


r 

x eq 


x 2 

d.o.f. 


Cu+Cu 
62.4 GeV 


0- 6% 
6 - 15% 
15 - 25% 
25 - 35% 
35 - 45% 


825 
681 
494 
340 
230 


3.2 
2.2 
1.4 
0.9 
0.4 


1.12 
1.11 
1.09 
1.09 
1.08 


3.70 
3.70 
3.70 
3.70 
3.70 


5.16 
4.92 
4.73 
4.57 
4.45 


4.7 
49 
2.4 
49 
1.5 
40 
3.0 
40 
4.9 
40 


Cu+Cu 
200 GeV 


0- 6% 
6 - 15% 
15 - 25% 
25 - 35% 
35 - 45% 
45 - 55% 


1474 
1129 
791 
536 
349 
216 


7.1 
6.2 
5.4 
4.9 
4.3 
4.2 


1.08 
1.07 
1.06 
1.05 
1.05 
1.04 


4.03 
4.03 
4.03 
4.03 
4.03 
4.03 


2.45 
2.40 
2.35 
2.31 
2.28 
2.24 


2.0 
49 

1.0 
49 
1.8 
49 

3.5 
49 

5.7 
49 
8.2 
49 


Au+Au 
19.6 GeV 


0- 6% 
6 - 15% 
15 - 25% 
25 - 35% 
35 - 45% 


1691 
1323 
966 
672 
429 




1.23 
1.19 
1.15 
1.12 
1.10 


2.90 
2.90 
2.90 
2.90 
2.90 




0.7 
28 
0.4 
28 
0.5 
28 
1.5 
28 
1.9 
27 


Au+Au 
130 GeV 


0- 6% 
6 - 15% 
15 - 25% 
25 - 35% 
35 - 45% 
45 - 55% 


4233 
3318 
2313 
1559 
1005 
615 


13.2 
11.9 
10.9 
10.0 
9.3 
8.6 


1.02 
1.02 
1.01 
1.01 
1.00 
1.00 


3.56 
3.56 
3.56 
3.56 
3.56 
3.56 


2.64 
2.45 
2.28 
2.14 
2.05 
1.93 


3.7 
49 
1.0 
49 
1.3 
41 
3.1 
41 
6.2 
41 
9.7 
41 


Au+Au 
200 GeV 


0- 6% 
6 - 15% 
15 - 25% 
25 - 35% 
35 - 45% 
45 - 55% 


5123 
3987 
2827 
1916 
1251 
762 


26.3 
24.8 
23.7 
22.7 
21.9 
21.1 


0.93 
0.93 
0.92 
0.92 
0.92 
0.91 


3.51 
3.51 
3.51 
3.51 
3.51 
3.51 


3.20 
3.08 
2.97 
2.90 
2.83 
2.76 


1.1 

49 

0.8 
49 

2.3 
49 
7.9 
49 

11.7 
49 

15.4 
49 
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Fig. 1: Dependence of the Diffusion- Model parameters for heavy systems (central Au + Au at RHIC 
energies) on the center-of-mass energy v / ijvjv: Quotient of interaction time and relaxation time 
for sinh- and exponential (dashed) extrapolation (upper frame); width of the peripheral sources 
including collective expansion (middle frame); effective width of the midrapidity source (lower 
frame). The results are for charged-hadron rapidity distributions, with extrapolations to LHC 
energies of 5.52 TeV. The dots refer to the fit values at ^sjvjv = 19.6, 130 and 200 GeV. 



is performed through the Jacobian 

J(y, = cosh(?7)- 

[l + ((m)/(p T )) 2 +sinh 2 ( ?? )]- 1 /2. (5 ) 

We approximate the average mass < m > of produced charged hadrons in the central region 
by the pion mass m^, and use a mean transverse momentum < pt > = 0.4 GeV/c. Due to 
the conversion, the partial distribution functions are different from Gaussians. 

In the linear two-sources version, the Relativistic Diffusion Model had been applied to 
pseudorapidity distributions of produced charged hadrons in Au+Au collisions at RHIC 
energies of 130 GeV and 200 GeV by Biyajima et al. [9]. However, it soon turned out from 
the net-proton results [6], and from general considerations, that an additional midrapidity 
source is required [6, 14]. This source for particle production arises mostly from gluon-gluon 
collisions and emerges at very short times. In the model we assume that it is generated at 
t=0 and y eq with full strength, and then spreads in rapidity space according to Eq.(3) during 
the strong-interaction time. It comes close to local thermal equilibrium with respect to the 
variable rapidity during the interaction time Tint and hence, we use the notion R l f q c (y, t) for 
the associated partial distribution function in y-space, with charged particles. 
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Fig. 2: Calculated pseudo-rapidity distributions of produced charged particles from Au + Au col- 
lisions at t/snn = 19.6, 130 and 200 GeV for six collision centralities (sequence as in Figure 3, 
bottom frames for central collisions) in comparison with PHOBOS data [15] . The analytical RDM- 
solutions are optimized in a fit to the data. The corresponding x 2 -values are given in Table 1. 
Dashed curves show the midrapidity sources for hadron production. 



The validity of this picture was underlined by our recent investigation of the d + Au 
system at 200 GeV in the three-sources model [10]. Here, an accurate modeling of the 
gradual approach of the system to thermal equilibrium in rapidity space was obtained. 
In particular, the dependence of the asymmetric pseudorapidity distribution functions on 
centrality was precisely described. In the present investigation, however, we concentrate on 
heavy symmetric systems where QGP formation is more likely. 

The dependence of the diffusion- model parameters on incident energy in central Au + Au 
collisions at RHIC is displayed in Figure 1 . Resulting values for the time parameter Tint /t v 
are shown as function of incident energy in the upper frame, with a functional dependence 
on the absolute value y ma x of the beam rapidity and hence, on energy given by 

• 1/ \ ( 6 ) 

T y smh(y max ) 

that is discussed in more detail in [17], whereas the dashed curve assumes an exponential 
dependence on log(y / S7v r /v)- The partial widths as functions of energy within the RHIC 
range for Au + Au arc shown in the middle and lower frames of Figure 1 for both peripheral 
and midrapidity sources, which differ for produced hadrons. The widths arc effective values: 
beyond the statistical widths that can be calculated from a dissipation-fluctuation theorem, 
they include the effect of collective expansion. Here we have plotted the values resulting 
from the ^-minimization that include the time evolution up to Tint, including collective 
expansion 

I^e, = [81n(2) • D\[{ eq ■ t v ■ (1 - exp(-2r mt /r,))] 1 / 2 . (7) 
The charged-particle distribution in rapidity space is obtained as incoherent superposi- 
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Fig. 3: Pseudorapidity distributions of charged particles from Cu + Cu collisions at ^/sjvjv = 62.4 
and 200 GeV for six collision centralities in comparison with PHOBOS data [16]. The analytical 
diffusion-model solutions are optimized in a fit to the data, with constraints discussed in the text. 
The corresponding ^-values are given in Table 1. 



tion of nonequilibrium and local equilibrium solutions of (3) 

dN ch (y,t = T int ) M i u t \ 

^ = N c h R l(V:^nt) 

+ N 2 ch R 2 ( Vl T int ) + N2R l e ° q C (y, T int ) (8) 

with the interaction time Tint (total integration time of the differential equation). In the 
present work, Ti n t/T y is determined together with the set of other free parameters Ti^eij and 
n ch f rom the x 2 — minimization with respect to the data and hence, the explicit value of Tint 
is not needed as an input. The resulting values for Ti n t/T y are given in Table 1. Although 
there are rather rapid changes of x 2 m narrow intervals of T in t / T y due to the simultaneous 
dependence on the other parameters, the fitting procedure provides a reliable criterion for 
the determination of the diffusion-model parameters. 

Heavy systems at RHIC. — Our results for pseudorapidity distributions of produced 
charged hadrons at six different centralities and three incident energies in Au + Au collisions 
are shown in Figure 2 in comparison with PHOBOS data [15]. Corresponding results for Cu 
+ Cu at two incident energies are given in Figure 3 compared with preliminary PHOBOS 
data [16] . At the lowest energy, only two sources are needed for the optimization of the RDM- 
parameters in a x 2 -fit, whereas three sources are indeed required at the higher energies. At 
the highest energy of 200 GeV, the Cu + Cu system requires a smaller percentage (7%) 
of particles in the midrapidity source compared to Au + Au, where 26% of the produced 
hadrons are in the equilibrated source for central collisions. This is consistent with the 
expectation that heavier systems are more likely to produce a locally equilibrated quark- 
gluon plasma. 
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The parameters of these calculations are summarized in Table 1 together with % 2 /d.o./.. 
The number of degrees of freedom (d.o.f.) is the number of data points minus the number 
of free parameters. We have not aimed at the absolute minimum of x 2 as in [9, 10] because 
this does not provide a sufficiently precise determination of the free parameters. Instead, 
we use physically meanigful constraints to reduce the number of free parameters, and to 
determine a local minimum. 

In particular, we impose a linear decrease of the time parameter T int /Ty, of the percentage 
of particles in the midrapidity source n^, and of the partial widths Pi, 2, eg with increasing 
impact parameter. Using these constraints, it is possible to obtain excellent results in 
the x 2 -optimization. Regarding the centrality dependence at fixed incident energy, the 
increase of the size of the midrapidity source towards more central collisions provides a 
good reproduction of the data. This is physically reasonable since the midrapidity source is 
expected to be more important towards more central collisions, where it may originate from 
an equilibrated parton plasma because of the high energy density. 

In a x 2 — minimization without any physical constraints [10], the results for the size of the 
central source as function of centrality have not shown such a trend. Instead the percentage 
of hadrons in the midrapidity source rises for more peripheral collisions, because the number 
of charged hadrons produced in the beam-like regions of pseudorapidity space falls more 
strongly than the overall number of produced particles. In view of the good quality of the 
constrained centrality dependence in the present work, however, this particular result which 
is physically difficult to understand may turn out to be an artifact of the fit procedure. 

Based on our present calculations at RHIC energies, we have also predicted rapidity 
distributions of produced charged hadrons in central Pb + Pb collisions at LHC energies 
using the extrapolations displayed in Figure 1. The results are shown in [17]. 

Conclusion. — To conclude, we have described charged-hadron pseudorapidity distri- 
butions in collisions of heavy symmetric systems at four RHIC energies and six centralitics 
with high precision in a three-sources Relativistic Diffusion Model (RDM) for multiparticle 
interactions. Analytical results for the pseudorapidity distribution of charged hadrons at all 
investigated centralities are found to be in excellent agreement with the available data. An 
extrapolation of the transport parameters to LHC energies has been performed. 

At the highest RHIC energy of 200 GeV, about 26% of the charged hadrons in central 
Au + Au collisions are produced in the midrapidity source. These particles come very close 
to statistical equilibrium in rapidity space during the strong-interaction time. In central 
collisions with very high energy density, they are likely to originate from a thcrmalized 
quark-gluon plasma. The midrapidity source is less pronounced towards more peripheral 
collisions, lower energies, and in smaller systems such as Cu + Cu. It vanishes at the lowest 
energy of 19.6 GeV. 
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